Dynamics of FRW Universes Sourced by Non-Linear Electrodynamics 
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We apply the dynamical systems tools to study the (linear) dynamics of Friedmann-Robertson- 
Walker universes that are fuelled by non-linear electrodynamics. We focus, mainly, in two particular 
models. In the first model the cosmic evolution is fuelled by cold dark matter, a cosmological 
constant and a non-linear electrodynamics field. In the second case non-singular cosmology and late- 
time accelerated expansion are unified in a model where the Einstein's field equations are sourced 
only by cold dark matter and a non- linear electrodynamics field. It is shown that, in contrast to 
previous claims, the cosmological effects coming from the non-linear electrodynamics field are not as 
generic as though. In fact, critical points in the phase space that could be associated with non-linear 
electrodynamic effects are not found. 

PACS numbers: 04.20.-q, 98.80.-k, 98.62. En, 98.80.Cq, 98.80.Jk 
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I. INTRODUCTION 



Studying the equations of the non-linear electrodynam- 
ics (NLED) is an attractive subject of research in general 
relativity (GR) thanks to the fact that such quantum 
phenomena as vacuum polarization can be implemented 
in a classical model through their impact on the prop- 
erties of the background space-time. Exact solutions of 
the Einstein's field equations coupled to NLED may hint 
at the relevance of the non-linear effects in strong grav- 
itational and magnetic fields. It has been speculated, in 
particular, that very strong electromagnetic fields might 
help avoiding the occurrence of space-time singularities 
[l[ . The impact of very strong electromagnetic fields (and 
of the NLED effects) regarding the causality issue in cos- 
mology is also of relevance @ . 

A different subject of research within the cosmologi- 
cal setting in GR coupled to NLED, is related with the 
chance for the NLED field to fuel primordial inflation. 
The cosmological inflationary scenario was proposed for 
the first time in the Ref. [|| where, in order to solve (or 
avoid) several problems of the standard model - such as 
the flatness and the horizon problems, amongts others -, 
it was anticipated that a self-interacting scalar field with 
a particular form of the self-interaction potential as the 
source of the Einstein's field equations, might cause the 
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universe to expand in an inflationary (super-accelerated) 
fashion. The inflationary paradigm has been supported 
by the observational evidence (ij. 

In the reference Q an anisotropic cosmological model 
coupled to Born-Infeld NLED was explored. It was found 
that this model might explain early time inflation Q. 
Thanks to the inflationary stage, the initially anisotropic 
universe might eventually isotropize. An alternative non- 
isotropic model where, additionally to the Born-Infeld 
NLED field, a cosmological constant was added to the 
Einstien's field equations, has been studied in Q ■ Yang- 
Mills cosmology with a non-Abelian Born-Infeld action 
has been also considered @]. 

Magnetic universes (vanishing electric component) 
have been investigated within the context of theories 
given by the Lagrangian density C\ = — \F + aF 2 + (3G 2 
(see below) Q . The non-linear term cx F 2 may cause the 
universe to bounce thus avoiding the initial (big-bang) 
sine ularity. Then, in Ref. [3 it was demonstrated that 
the inclusion of a non-vanishing electric component E 
where E 2 ~ B 2 , removes the bounce and, in consequence, 
the universe starts its evolution in a singular state. Mod- 
els with Lagrangian density of the form £2 = ~jF — 
may account for the late-time stage of accelerated ex- 
pansion of magnetic universes [ll[ . In the later case, if a 
non-vanishing electric component E is considered, accel- 
erated expansion is allowed only when E 2 < 3B 2 [nj. 

Aim of the present investigation is to explore the 
asymptotic properties of cosmological models where gen- 
eral relativity is coupled to NLED given by the above 
Lagrangian densities C\ and C%. We will rely on the 
use of the dynamical systems tools with the hope that 
such relevant concepts as past and/or future attractors, 
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or saddle equilibrium points, could be correlated with 
generic cosmological behaviour. Our goal will be to write 
the cosmological (Einstein's) equations in the form of 
an autonomous system of ordinary differential equations 
(ODE). Although these equations are non-linear, we will 
expand them in the neighbourhood of the equilibrium 
(critical) points up to the linear approximation. Eval- 
uating the sign of the real parts of the eigenvalues of 
the corresponding linearization matrices will allow us to 
judge about the stability of the critical points. We will 
show, in particular, that equilibrium points correspond- 
ing to dominance of NLED effects are not found, signal- 
ing at a fundamental drawback of the models. Here we 
use natural units (8irG = Sir /m 2 Pl = H = c = 1). 

II. A-CDM-NLED MODEL 

In this and in the next sections we shall focus in a FRW 
cosmological model with flat spatial sections, fuelled by 
three generic sources: i) pressureless cold dark matter 
(CDM), ii) dark energy in the form of a non- negative cos- 
mological constant A > 0, and iii) non- linear electrody- 
namics (NLED). While the first and second components 
(CDM and A-dark energy, respectively) are required for 
consistency of the late-time cosmological evolution, the 
third component (properly the non-linear electrodynam- 
ics field) is important at early as well as at intermediate 
times. At early times departure from standard radiation- 
like behaviour is apparent, however, later on in the course 
of the cosmic evolution the NLED field behaves just like 
a standard radiation field. 

In the next section we want to explore the asymptotic 
properties of the above model that could be correlated 
with generic cosmological behaviour. This will allow us 
to make a judgment about the actual relevance of NLED 
on early-time dynamics, in particular for the bounce of 
the universe. Later on, in section IV, we shall explore a 
model where the bounce and late-time inflationary dy- 
namics are originated by a single NLED field. 

In general the Lagrangian density of NLED is given 
by C = C(F, G) where F = F^F^ and G = F^F^ 
are the electromagnetic invariants (F^u = \^^vaf)F al3 ^ . 
The energy-momentum tensor that is associated with this 
Lagrangian can be written in the following way: 

T„„ = -l{d F L)F«F av + (Gd G C - C)g^. (1) 

In order to meet the requirements of homogeneous and 
isotropic cosmologies (as, in particular, the one associ- 
ated with FRW space-times), the energy density and the 
pressure of the NLED field should be evaluated by av- 
eraging over volume (for details see [12j and references 
therein). Additionally it has to be assumed that the 



electric and magnetic fields, being random fields, have 
coherent lengths that are much shorter than the cosmo- 
logical horizon scales. After several considerations the 
energy-momentum tensor of the electro-magnetic (EM) 
field associated with the Lagrangian density C — C(F, G) 
can be written in the form of the energy-momentum ten- 
sor for a perfect fluid: 

T^y = (p+p) Ufj,u v - pg^, (2) 

where 

p = -£ + Gd G C - 4d F CE 2 , (3) 

P = C- Gd G C - ± (2B 2 - E 2 ) d F C, (4) 

E and B being the averaged electric and magnetic fields 
respectively. For the purposes of the present investigation 
we shall consider the following Lagrangian density that 
has been formerly studied in |9j: 

£=-^F + aF 2 + pG 2 , (5) 

where a and f3 are arbitrary constants. As mentioned 
in [l2j], a particularly interesting case arises when only 
the average of the magnetic part B is different from zero, 
leading to the so called magnetic universe. This case 
turns out to be relevant in cosmology as long as the av- 
eraged electric field E is screened by the charged pri- 
mordial plasma, while the magnetic field lines are frozen 
[ijl] . Let us focus here, exclusively, in a non-singular uni- 
verse proposed in Ref.[l3| (Lagrangian density (O with 
/3 = 0). According to this scenario the energy density ([3]) 
and effective presure Q are given by: 

PB = \b 2 (I - 8aB 2 ) , (6) 

PB = \b 2 (1 - 40aB 2 ) , (7) 

respectively, where, in order to produce the bounce, the 
magnetic field B has to be associated with the following 
evolution law: B — (Bo is a constant of integration) . 

Notice that the energy density of the NLED field B 
vanishes at the maximum allowed value of the magnetic 
field B = l/y/8a. ps is a minimum at B — 0, while at 
B = l/\/16a it is a maximum instead. At the maximum 
the following relationship between the energy density and 
the parametric pressure arises: pb = ~Pb, i- e., at the 
maximum the NLED field behaves as vacuum does. In 
consequence the equation of state parameter lob fluctu- 
ates between 1/3 (radiation) and —1 (vacuum energy). 



3 



This means, in turn, that the NLED field could con- 
tribute towards primordial (early-time) inflation. How- 
ever, it will be shown here that this inflationary be- 
haviour is not as generic as though since it can not be 
correlated with any equilibrium point in the phase space 
of the model. Notice that the energy density is positive 
(ps > 0) whenever 
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B < 

2V2a' 

meanwhile the pressure will be negative only if 
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B > 



The strong energy condition is violated (ps + 3ps < 0) 
whenever 
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B > 

In the later case the NLED field might fuel early-time 
inflationary dynamics. As already said, it will be demon- 
strated that inflationary NLED does not represent a crit- 
ical point in the phase space of the model, meaning that 
it is unlikely that a NLED field might be a generic source 
of cosmic inflation. 

The starting point for the present investigation will be 
the following Einstein's field equations (as usual H = - 
is the Hubble parameter, the overdot accounts for deriva- 
tive in respect to the cosmic time t, while ps and ps are 
given by © and ([7]) respectively): 



3H 2 = pc dm + A + p B , 
2i? = -peon - (pb + Pb) , 

plus the following conservation equations: 



(8) 



PCDM + SHpcDM = 0, (9) 

and the continuity equation for the NLED magnetic field, 



Pb+3H( Pb +Pb) = 0. 



(10) 



Equations (J8J) , (|9J) and (|10| form the mathematical basis 
of the model of interest in the present paper. 



III. DYNAMICAL SYSTEMS STUDY 

The dynamical systems tools offer a very useful ap- 
proach to the study of the asymptotic properties of the 
cosmological models [3]. In order to be able to apply 
these tools one has to (unavoidably) follow the steps enu- 
merated below. 

• First: to identify the phase space variables that 
allow writing the system of cosmological equations 
in the form of an autonomus system of ordinary 
differential equations (ODE). There can be several 
different possible choices, however, not all of them 
allow for the minimum possible dimensionality of 
the phase space. 

• Next: with the help of the chosen phase space vari- 
ables, to build an autonomous system of ODE out 
of the original system of cosmological equations. 

• Finally (some times a forgotten or under- 
appreciated step): to indentify the phase space 
spanned by the chosen variables, that is relevant 
to the cosmological model under study. 

After this one is ready to apply the standard tools of 
the (linear) dynamical systems analysis. 



A. The Autonomous System of ODE 

We introduce the following phase variables: 



B 



.... Va 

V V3H' 



11) 



We will be focused, exclusively, on expanding FRW uni- 
verses, so that x > 0, y > 0, and v > {B > 0). 

After the above choice of variables of the phase space, 
the system of cosmological equations (I8I9I10[) can be 
translated into the following autonomous system of or- 
dinary differential equations (ODE): 



H' 



V 



-II 



W 1 



IT' 



f = -^[l-2/ 2 + y(l-15xV)], (12) 

where the prime denotes derivative with respect to the 
new time variable r = In a - properly the number of efold- 
ings. The phase space where to look for the equilibrium 
points of the above system of ODE is the following: 
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* = {(x, y, v) : < y 2 + x 2 (l - 3x 2 v 2 ) < 1, 

x > 0, y > 0, v > 0, x 2 v 2 < 1/3}. (13) 

Additionaly it will be helpful to have the following 
parameters of observational importance Qb = Pb/3H 2 
- the NLED field dimensionless energy density parame- 
ter -, and the equation of state (EOS) parameter = 
Pb/ PBi written in terms of the phase space variables: 

n B = x 2 (l - 3arV), (14) 

1 f l-lbx 2 v 2 \ 
" B = 3 U-3sV J' (15) 



fl A = y 2 , Qcdm = 1 - A - Rb. (16) 

Notice also, that for the deceleration parameter one 
has q = -(l + H'/H). 

B. Critical Points 

The critical points of the autonomous system of equa- 
tions (JT2J) and the eigenvalues of the corresponding lin- 
earization (Jacobian) matrices, are summarized in table 
III 

There are three equilibrium points at all: i) CDM- 
dominated solution (point Pi in TabHJ, ii) radiation- 
dominated solution (point P2), and iii) de Sitter - in- 
flationary - solution (point P3). The latter is a non- 
hyperbolic critical point so that we are not able to make 
final jugment about its stability. Maximum we can say 
is that there is a stable subspace attached to P 3 , that is 
spanned by eigenvectors associated with Ai and A2. How- 
ever, the plots of the orbits of the autonomous system of 
ODE (fT2|) shown in the figure Fig[IJ demonstrate that the 
de Sitter solution is the late-time (future) attractor in the 
phase space of th e mod el. The fact that, at the attrac- 
tor point v — Ay/ aA/3 — > 0, means that the non-linear 
electrodynamic effects - represented by the constant pa- 
rameter a - have to be very weak. The remaining points 
Pi and P2 are always saddle critical points in ^. There is 
not any past attractor for the orbits of the system of ODE 
(TTZ|) . From the cosmological point of view the model is 
viable since there are critical points associated with ra- 
diation and CDM dominance, as well as with the present 
stage of accelerated expansion. The same results would 
arise if one replaced the non-linear electrodynamics field 
by a standard Maxwell one. 




FIG. 1: The orbits of the autonomous system of ODE (fl2jl 
for arbitraryly chosen sets of initial conditions (upper panel) . 
Only expanding cosmologies are reflected in the phase space. 
The projection of the corresponding orbits in the phase plane 
(x, y) are shown in the left hand lower panel, while the flux 
in time of the orbits corresponding to expansion are shown 
in the right hand lower panel. Notice that the de Sitter solu- 
tion (0, 1, 4v/ aA/3) is the late-time (future) attractor of these 
orbits, while the solution dominated by the CDM (0, 0, 0) is 
always a saddle critical point in the phase space. 



It comes as a surprise that the non-linear electrody- 
namics effects - representing the strong field limit where 
the scalar of curvature is small and the volume of the uni- 
verse attains its minimum [l5| -, can not be associated 
with any equilibrium point in the phase space ^ (|13[) . 
This means that the NLED effects are not as generic 
as previously though. In particular, the non-singular be- 
haviour correlated with the bounce, is not a generic prop- 
erty of the model where the radiation field is replaced by 
a (magnetic) NLED field. The bounce is just a particular 
solution of the model that, depending on the initial data 
chosen, might or might not occur, i. e., it could be only 
a unstable state in the phase space. 

In the next section we shall explore a model where also 
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TABLE I: Properties of the critical points for the autonomous system (I12|l . and eigenvalues of the linearization matrices. 
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&CDM 
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Pi 











1 








1/3 


1/2 


-3/2 


-1/2 


3/2 


Pi 


1 
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1/3 
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-2 


1 


2 


Pz 
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4,/aA/3 





1 





1/3 


-1 


-3 


-2 






the fate of the cosmic expansion is determined by a NLED 
field (no A-term at all). We shall show, in particular, 
that from the cosmological point of view this model is 
not attractive, since, in this case there are no critical 
points in the corresponding phase space that could be 
identified with late-time cosmological dynamics. 

IV. UNIFIED DESCRIPTION OF BOUNCING 
COSMOLOGY FOLLOWED BY LATE-TIME 
ACCELERATED EXPANSION 

In this section we focus in the study of the asymp- 
totic properties of a cosmological model with interesting 
features, namely a phase of current cosmic acceleration 
and the absence of an initial singularity, which was pro- 
posed in [13] (see [ll|) and is based upon the following 
Lagrangian density: 

C = -^F + aF 2 +7F- 1 , (17) 

where 7 is a new overall (constant) parameter. The sec- 
ond (quadratic) term dominates during very early epochs 
of the cosmic dynamics, while the Maxwell term (first 
term above) dominates in the radiation era. The last 
term in (fl7| is responsible for the accelerated phase of 
the cosmic evolution [12], so that here we can safely ig- 
nore the cosmological constant term A in Einstein's field 
equations |(5J). The above Lagrangian density yields a 
unified scenario to describe both the acceleration of the 
universe (for weak fields) and the avoidance of the ini- 
tial singularity, as a consequence of its properties in the 
strong-field regime. 

After averaging over a (flat) FRW space-time the 
stress-energy tensor associated with (flTf can be written 
in the form of an equivalent perfect fluid stress-energy 
tensor with energy density and parametric pressure: 



respectively. Notice that, for large values of the NLED 
field, positivity of energy requires that B < l/\/8a, 
while, for small enough values of B <C 1, positivity of 
energy implies that B > 7 1 / 4 . Therefore, in the uni- 
fied model of [12] (Lagrangian density (|17p) there arise 
both higher and lower bounds on the values the NLED 
field B can take. However, given that the observa- 
tional data constraints the parameter 7 to be v/pyf ~ 
4 x 10~ 28 g cm~ 3 [ll[, then the lower bound on B can be 
admitted without going into conflicts with observations. 

Our goal in this section will be to put the cosmolog- 
ical equations {l}, (|9|), and (fl"0|) . with ps and ps given 
by (fl~8"|) , in the form of an autonomous system of ODE. 
Recall that, since the NLED Lagrangian density (jTTJ) ac- 
counts also for the present accelerated stage of the cosmic 
expansion, then in the Einstein's field equation ([5]) the 
cosmological constant is set to zero A = 0. In the present 
case we choose the following phase space variables: 



X =-jte> " = (19) 

The corresponding autonomous system of ODE can be 
written in the following form: 

x =-x\2+ — I , 

V a) 

H> „ H' 



H' _ 3 fx i -15x 6 v 2 + 7w 2 
IT ~ ~2 ~ V 2^ 




The phase space relevant for this case can be defined as: 

$>u = {(x, v, w) : < x 4 - 3x 6 v 2 - w 2 < x 2 , 

x>0,v>0,w>0}. (21) 



PB = ^(l-8aB 2 -^\ 



2 \ B 4 J ' 

PB = |- (l - 40aB 2 + |l) , (18) 



In terms of the above phase space variables: 
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TABLE II: Properties of the critical points for the autonomous system (|20[) . and eigenvalues of the linearization matrices. 
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&CDM 
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Ai 


A 2 


A 3 


Pi 











1 





1/3 


1/2 


-3/2 


-1/2 


3 


p> 


1 











i 


1/3 


1 


-2 


1 


4 



.x — 3x v 



6, ,2 



x 4 — I5x 6 v 2 



7w' 



3x 6 v 2 



(22) 



There can be found only two critical points of the au- 
tonomous system of ODE ([20]) - shown in Tab [III i) the 
CDM-dominated solution (point Pi in Tab HH) , and ii) 
the radiation-dominated solution (critical point P2). As 
it is apparent from Tab [III both represent saddle equi- 
librium points in the phase space (|2 1 1) . Notice that 
there is neither past nor future attractors that could be 
associated with the model of this section. This fact is 
specially illustrated in the figure FigO where the orbits 
of the system of ODE ([20)) are plotted in a portion of the 
3D phase space (x,v,w) S ^u- 

In the present model, as in the model of section III, the 
possibility that the NLED field might cause the late-time 
accelerated expansion to happen - although existing -, 
is not a generic property of the model, since there are 
no critical points that could be correlated neither with 
early-time nor with late-time NLED effects. This fact 
is reflected in the figure Fig. 2, where it is apparent that 
the orbits of the system of ODE evolve from points that 
are dominated by early-time non-linear effects (large v-s 
and vanishing w-s) towards larger values of the variable 
w, that amounts to increasing relevance of the last term 
in the NLED Lagrangian density C of (JT7J) . However, as 
clearly seen, there is neither any past (attractor) source 
points from which the orbits are repelled in the past, nor 
converging focus towards which the orbits are attracted 
in the future. 




FIG. 2: The orbits of the autonomous system of ODE ([20)1 
for arbitraryly chosen sets of initial conditions (upper panel) . 
Only expanding cosmologies are reflected in the phase space. 
Different perspectives of the flux in time of the orbits are 
shown in the lower panels. Notice that the solution dominated 
by the CDM (point (0, 0, 0) in Tab HI)) is always a saddle point. 



V. DISCUSSION AND CONCLUSIONS 

Non-linear electrodynamics can supply a useful sce- 
nario where to discuss two relevant problems of the stan- 
dard cosmological model: i) the initial singularity related 
with the big-bang, and ii) the present stage of accelerated 
expansion of the universe. A Lagrangian density of the 
form given in ([5]) might provide a solution to the first 
problem, since in this model the universe undergoes a 



bounce without singularity of the curvature invariants. 
There can be also a chance for this model to fuel the 
early-time inflation. A model based on the Lagrangian 
(fTT)) . instead, can account for the present accelerated 
stage of the cosmic evolution, as well as for avoidance 
of the initial singularity, in a unified picture. 

A strightforward analysis of the ordinary differential 
equations ()12j) or ([2"0")) reveals, however, that there are no 
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critical points in the phase space of the above models that 
could be associated with non-linear effects. Actually, in 
both models, existence of equilibrium points correlated 
with non-linear effects requires that x 1 = v' = w' = 
and, simultaneously x 7^ 0, v ^ 0, and w ^ 0. 1 However, 
a simple inspection of the corresponding equations in (|12|) 
and (j20|) reveals that the simultaneous fulfillment of the 
above conditions is imposible. 

The knowledge of the equilibrium points in the phase 
space corresponding to a given cosmological model is a 
very important information since, independent on the ini- 
tial conditions chosen, the orbits of the corresponding 
autonomous system of ODE will always evolve for some 
time in the neighbourhood of these points. Besides, if 
the point were a stable attractor, independent on the 
initial conditions, the orbits will always be attracted to- 
wards this point (either into the past or into the future). 
Going back to the original cosmological model, the ex- 
istence of the equilibrium points can be correlated with 
generic cosmological solutions that might really deside 
the fate and/or the origin of the cosmic evolution. In 
a sense the knowledge of the asymptotic properties of 
a given cosmological model is more relevant than the 
knowledge of a particular analytic solution of the cor- 
responding cosmological equations. While in the later 
case one might evolve the model from given initial data 
giving a concrete picture that can be tested against exist- 
ing observational data, the knowledge of the asymptotic 
properties of the model gives the power to realize which 
will be the generic behaviour of the model without solv- 
ing the Einstein's field equations. In the dynamical sys- 
tems language, for instance, a given particular solution 
of the Einstein's equations is just a single point in the 
phase space. Hence, phase space orbits show the way the 
model drives the cosmological evolution from one partic- 
ular solution into another one. Equilibrium points in the 
phase space will correspond to solutions of the cosmolog- 
ical (Einstein's) equations that, in a sense, are preferred 
by the model, i. e., are generic. The lack of equilibrium 
points that could be correlated with a given analytic so- 
lution of the model, amounts to say that this solution is 
not quite generic, otherwise unstable in terms of phase 
space variables, and can not be taken too seriously. 

This is precisely the case with the models under study 
in the present paper. The analysis of the Einstein's 
field equations, together with the conservation equations, 
shows that, in principle, it is possible that the effects of 
NLED fields might originate the early-time inflation as 
well as the late-time stage of accelerated expansion of the 
universe. The possibility to avoid the initial singularity, 



replacing it by a bounce, seems to be another relevant 
property of the model. However if one applies the tools 
of the dynamical systems approach, one can realize that 
this is no more than just a hope. Actually, solutions 
where the big-bang singularity is replaced by a bounce, or 
where the NLED effects fuel primordial inflation and/or 
late-time accelerated dynamics are not actually preferred 
by the model. Instead, the solutions dominated by CDM 
and by standard radiation represent equilibrium points 
in the phase space, meaning that the cosmic evolution 
might transit through the corresponding stages during a 
quite long time. A clear deficiency of the models is the 
inexistence of past attractors, meaning that the initial 
conditions from which the universe might be evolved (in 
the vicinity of the bounce) are not generic neither. The 
inexistence of late-time attractors in the second model 
(unified description of bouncing cosmology with late-time 
accelerated expansion), means that the fate of the cos- 
mic evolution in the model, although inflationary, is quite 
uncertain. 

Summarizing we can say that the second model, where 
the bounce and the late-time accelerated expansion are 
described in a unified picture, although at first sight 
seems to be a good candidate to address the initial sin- 
gularity and the late-time inflation problems at once, in 
fact can not be a competitive model since both phenom- 
ena are not as generic as expected, so that it suffers from 
such serious drawbacks as the fine tunning and the co- 
incidence problems. The A-CDM-NLED model, instead, 
is a good candidate to address the late-time acceleration 
of the expansion thanks to the A-term. In general the 
cosmic evolution in this model transits by the radiation- 
dominated followed by the CDM-dominated (intermedi- 
ated) stages - corresponding to equilibrium points in the 
phase space of the model -, to eventually evolve into the 
de Sitter solution 3H 2 = A. However, the behaviour that 
is determined by the NLED term in this case - the bounce 
- can not be associated with any equilibrium point in the 
phase space, meaning that, as in the second model, the 
NLED effects are not generic and can be met just as a 
particular solution of the Einstein's field equations. 
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1 The variable v accounts for the influence at early timesof the by 7 in (Q2} which dominates at late times. 

NLED term that is multiplied by the parameter a in JT7} - as 
well as in J5j — , while w is correlated with the term multiplied 
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